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ABSTRACT 

A  comparison  is  made  of  different  coherent  optical 
filters  with  respect  to  their  ability  to  locate  an  edge  and 
identify  phase  information.  These  different  filtering  methods 
are  bandpass  filtering,  the  Schlieren  method,  the  Hilbert 
transform  method,  a  linear  amplitude  filter,  and  a  differential 
filter;  these  filters  are  one -dimensional.  The  differential 
filter  is  a  combination  of  two  filters,  a  linear  amplitude 
filter  and  a  Hilbert  transform  filter.  A  linear  amplitude 
filter  is  made  photographically  by  a  programmed  exposure  of 
film  to  give  the  correct  density  variations.  A  Hilbert  transform 
filter  is  a  phase  filter  which  is  accomplished  with  a  dielectric 
coating  on  glass.  Bandpass  filtering  is  accomplished  by  using  a 
wire.  With  the  use  of  a  razor  blade,  the  Schlieren  method  is 
abhieved.  The  comparison  of  these  filters  indicates  that 
bandpass  filtering  is  superior  for  the  location  of  an  edge  and 
that  the  Hilbert  transform  method  is  best  for  identifying  phase 
information. 
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max 

A  . 

.  min 

A(a)) 
A'  (oj) 
b 
B 


fcs 
F(aj) 
F'  ({») 
SF 


NOMENCLATURE 

Lower  cutoff  frequency  of  the  bandpass  filter 
Maximiim  value  of  the  filter  A'  (w) 

Minimum  value  of  the  filter  A' (u) 

Linear  amplitude  filter 
Approximate  linear  amplitude  filter 
Upper  cutoff  frequency  of  the  bandpass  filter 
Lower  cutoff  frequency  of  the  filter  A* (w) 

Lower  magnification  cutoff  frequency  of  the  filter 
A'  (oj) 

Lower  radial  cutoff  frequency  of  the  filter  A'  (to) 

Velocity  of  a  light  wave  in  a  vacuum 

_2 

Centimeters  (10  meters) 

Half  of  the  width  of  the  filter  A* (u) 

Film  density 

Normal  thickness 

Thickness  at  an  angle  of  45° 

Base  of  the  natural  logarithm 
Exposure 

Exposure  for  the  density  equal  to  zero 
Focal  length  of  a  lens 
Foot  candle  •  seconds 
Differential  filter 
Approximate  differential  filter 
Indicates  Fourier  transform 


vii 


AEDC-TR-71-137 


G(a)j^/Wy) 

G(x) 

g{x,y) 

H(io) 

I 


I. 

,  1 


I 


o 


j 


Indicates  inverse  Fourier  tremsform 
Fourier  transform  of  g(x,y) 

Amplitude  of  edge  due  to  bandpass  filtering 
General  variable 
Hilbert  transform  filter 
Intensity 

Intensity  in  the  image  plane 
Constant  intensity 


k 

k' 


k 


1 


L 

m 

mm 

mv 

MHz 

n 

”a 

"b 

nsec 

q 

RPM 

rpm 


Slope  of  the  filter  A'  (to) 

Slope  of  the  filter  A* (x^) 

Slope  of  the  linear  portion  of  the  amplitude 
transmission  versus  exposure  curve 
Slope  of  the  filter  F* (w^) 

Lens 

Magnification 

_3 

Millimeters  (10  meters) 

Millivolts  (lO”^  volts) 

Megahertz  (10^  hertz) 

Index  of  refraction 
Index  of  refraction  of  air 

Index  of  refraction  of  a  dielectric  material 
-9 

Nanoseconds  (10  seconds) 

Image  distance  from  lens 

Speed  of  the  D.  C.  motor  with  units  rpm 

Revolution  per  minute 
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Si 

t 

t' 


T(uj) 

U(x^) 

U(a)) 

V 


V 

o 


vdc 
V, 


DC 


W 


W_ 


m 


Amplitude  transmission  of  the  step  function 

Sine  integral 

Time 

Time  required  for  a  light  wave  to  travel  through 
a  certain  dielectric  material  thickness 
Amplitude  transmission 

Amplitude  transmission  for  the  exposure  equal  to 
2ero 

Filter  for  the  Schlieren  method 
Amplitude  in  the  image  plane 
Amplitude  in  the  frequency  plane 
Velocity  of  the  edge 
Slope  of  the  function  x^^^ 

Slope  of  the  function  x^2 

The  value  at  which  the  slope  of  x^^  and  x^2 
equal 

D.  C.  volts 

The  D.  C.  voltage  input  to  the  motor 
Upper  cutoff  frequency  of  the  filter  A' (w) 

Upper  magnification  cutoff  frequency  of  the 
filter  A'  (oj) 

Upper  radial  cutoff  frequency  of  the  filter  A'  (w) 
Distance  in  the  frequency  plane 
Distance  variable  for  the  portion  of  the  filter 
A' (x^)  described  by  the  H  and  D  curve 
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x^2  Distance  variable  for  the  portion  of  the  filter 

A' (x^)  described  by  the  amplitude  transmission 
versus  exposure  curve 
x^  x-coordinate  in  the  image  plane 

x^  x-coordinate  in  the  object  plane 

y-coordinate  in  the  image  plane 
y^  y-coordinate  in  the  object  plane 

y  Slope  of  the  linear  portion  of  the  H  and  D  curve 

6  (oj)  Delta  function 

C  Constant  angle  of  the  filter  F' (w) 

n  Positive  integer 

6  Angle  between  the  x  axis  of  a  filter  and  axis 

of  the  frequency  plane 
0(co)  Phase  filter 

ic  Wave  number 

X  Wavelength  of  light  source 

ysec  Microseconds  (10  ^  seconds) 

V  Velocity  of  a  light  wave  in  a  medium 

Velocity  of  a  light  wave  in  air 
Vj^  Velocity  of  a  light  wave  in  a  dielectric  material 

T  Width  of  the  function  s' (x^) 
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*  Convolution  integral 

''  Hilbert  transform 
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CHAPTER  I 

INTRODUCTION 

A  major  product  of  scientific  and  engineering 
endeavors  is  the  data  resulting  from  experiments  and  tests. 
In  this  thesis,  an  area  of  data  processing  which  might 
appropriately  be  called  data  enhancement  is  addressed.  Much 
of  the  data  encountered  in  this  work  has  arisen  from  aero¬ 
dynamic  simulation  and  testing;  however,  the  techniques 
which  will  be  discussed  are  not  limited  to  problems  in  the 
field  of  aerodynamics. 

Of  primary  interest  here  are  data  relating  to 
dimensional  changes  and  phase  changes.  Dimensional  changes 
might  correspond  to  the  ablation  of  a  projectile  fired 
through  a  chamber,  where  the  phase  changes  might  correspond 
to  density,  temperature  and  pressure  gradients  about  a 
subject  under  test.  The  techniques  which  are  employed  in 
this  investigation  are  coherent  optical  spatial  filtering 
techniques.  Basic  to  the  application  of  such  techniques  is 
the  simplicity  resulting  from  use  of  a  simple  positive  lens 
to  effect  a  two-dimensional  Fourier  transformation.  An 
imaging-transform  configuration  is  depicted  by  Figure  1 
which  may  be  employed  for  performing  spatial  filtering 
operations.  Such  a  system  utilizes  the  laser  as  its  source 
of  coherent  light. 
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Figure  1.  A  single  coherent  optical  data  processing  system 
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As  stated  above,  coherent  optical  data  processing 
Is  based  on  the  Fourier  transform.  The  Fourier  transform  Is 
defined  as 


®{g(x,y)}  = 


00  00 


J  J  g(x/ 


y)e 


-j(w  X  +  w  v) 


^  dxdy  =  G(u  ,u)  ) 
X  y 


(1) 


—  CO  —00 


and  the  Inverse  Fourier  transform  Is 


2F“^{G(w  )}  = 

X  y 


OO  GO 


4Tr 


iff  +j  (“xX  +  U)  y) 

^  J  J  G(w„/W„)e  *  y  dw„dw. 


X'  y' 


X  y 


(2) 


—  CO  —00 


A  positive  lens  will  perform  the  Fourier  transform  of  an 

Input  In  the  object  plane  of  the  lens.  The  transform  Is 

/ 

located  In  the  focal  plane  or  generally  referred  to  as  the 
frequency  plane.  The  frequency,  w.  Is  a  spatial  frequency 
of  units  llnes/mm.  The  Image  plane  of  the  lens  Is  the 
Fourier  transform  of  the  frequency  plane,  or 

SF  {3F{g(x^,y^) }}  =  g(x^,y^)  =  gC-x^.-y^)  (3) 

where  (x^/y^^)  are  the  coordinates  of  the  Image  plane,  and 
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(x^,y^)  are  the  coordinates  of  the  object  plane.  Both 

coordinate  units  are  usually  millimeters.  Generally,  the 

processing  of  data  in  the  object  plane  is  accomplished  by 

filtering  in  the  frequency  plane,  and  the  result  is  observed 

in  the  image  plane  [1].^ 

In  the  past,  a  variety  of  different  means  has 

been  used  to  process  optical  information  that  is  not  visible 

to  the  unaided  eye,  and  to  minimize  the  error  in  the 

location  of  an  edge  of  an  object.  Some  of  the  common 

methods  previously  used  are  bandpass  filtering,  the 

Schlieren  method,  the  Hilbert  transform  method,  and  the 

differential  filtering  method. 

The  Schlieren  method  is  effected  by  blocking  half 

* 

of  the  frequencies  in  the  frequency  plane  with  a  razor 
blade,  for  example.  The  Schlieren  method  may  be  described 
as  a  filtering  process,  the  filter  being  expressed  as  [1] 

T(u)  =  j(l  +  sgn(u)))  (4) 

where  the  function  sgn(u)  is  defined  as 

-1  oj  <  0 

sgn(u))  =  0  w  =  0  (5) 

1  0)  >  0 


^Numbers  in  brackets  correspond  to  similarly 
numbered  references  in  the  bibliography. 
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The  Schlieren  method  is  used  in  both  accentuating  an  edge  of 
an  input  and  in  the  conversion  of  phase  information  to 
intensity  information. 

The  Hilbert  transform  method  consists  of  applying 
to  the  negative  frequencies  a  phase  shift  of  +  ■rr/2  radians 
and  to  the  positive  frequencies  a  phase  shift  of  -  tt/2 
radians.  The  Hilbert  transform  filter  is  described  by  [2] 

H(£»))  =  -j  sgn(w)  (6) 

2 

where  j  =  -1.  This  filter  is  fabricated  by  coating  one- 
half  of  a  glass  plate  with  a  dielectric  material  to  give  a 
phase  shift  of  tt  radians.  The  thickness  of  the  dielectric 
material  retards  the  light  wave  by  X/2.  The  Hilbert  trans¬ 
form  is  used  to  convert  phase  information  to  intensity 
information  [3] . 

A  bandpass  filter  is  the  sum  of  a  high  pass  filter 
and  a  low  pass  filter.  The  high  pass  filter  is  realized  by 
placing  an  opaque  disc  -at  the  center  of  the  frequency  plane 
for  a  two-dimensional  input  or  a  wire  at  the  center  of  the 
frequency  plane  for  the  one-dimensional  input.  The  low  pass 
filter  may  also  be  realized  by  limiting  the  aperture  to  cut 
off  the  higher  frequencies.  In  practice,  the  low  pass 
filter  may  result  from  the  lens  used  to  perform  the  Fourier 
transform  of  the  input  in  the  object  plane.  The  bandpass 
filter  is  used  to  accentuate  the  edge  of  an  input. 
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The  differential  filter  is  a  complex  linear  ampli¬ 
tude  filter  of 


F(w)  =  jtJ 

=  |u|  j  sgn(w)  (7) 

The  differential  filter  is  composed  of  two  filters,  a  linear 
amplitude  filter  and  a  phase  filter.  This  theoretical 
filter  cannot  be  realized  because  it  is  physically 
impossible  to  have  F(u))  =  0  at  w  =  0  [4].  However,  an 
approximate  one-dimensional  differential  filter  can “be 
fabricated.  The  reasons  for  the  approximation  and  the  limi¬ 
tations  that  they  impose  are  discussed  in  Chapter  II.  Both 
phase  information  and  edges  of  an  input  are  accentuated  by 
the  differential  filter. 

Use  of  a  differential  filter  was  suggested  by 
Eguchi  and  Carlson  [5] .  Other  authors  have  stated  that  they 
used  a  differential  filter  in  optical  correlation  systems  to 
improve  the  signal-to-noise  ratio,  but  no  mention  was  made 
of  the  approximation  made  to  represent  a  differential  filter 
[6,7].  Some  authors  have  synthesized  optical  information 
and  a  differential  filter  on  the  computer  [8,9]. 

The  apparatus  and  the  methods  used  in  developing 
the  linear  amplitude  filter  and  phase  filter  are  discussed 
in  Chapter  II.  Chapter  III  gives  an  experimental  comparison 
of  the  various  optical  filtering  methods  operating  on  a 
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one-dimensional  step  function  and  phase  information. 
Comparative  conclusions  relevant  to  the  advantages  of  the 
various  filters  are  made  in  Chapter  IV. 
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CHAPTER  II 

THE  OPTICAL  DIFFERENTIAL  FILTER 

I.  THE  THEORY  OF  THE  DIFFERENTIAL  FILTER 

Consider  the  one-dimensional  case  of  an  input, 
g (x) ,  in  the  object  plane  of  the  lens  in  Figure  1,  page  2. 
The  Fourier  transform  of  g(x)  is  G(w).  To  find  the  type  of 
filter  needed  to  perform  the  derivative  of  the  input,  g(x) , 
the  Fourier  transform  of  the  derivative  must  be  known.  The 
Fourier  treuisform  of  dg/dx  is  [2] 


_  a 

dx  3x7 

o  o 


{G(u)}] 


_  1  d 

~  Tn  dx_ 
o 


f 


G(uj)  e 


Dwx, 


da>] 


jwx 

ja)G(a))e  du 


Therefore, 


(8) 


dg  (X  ) 

^  °  >  =  jwGCoj)  (9) 

o 


Thus,  in  order  to  differentiate  the  input,  a  filter  of 
complex  amplitude  transmittance,  ju,  is  placed  in  the 
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frequency  plane.  Then  the  function  in  the  Image  plane  is 

dg(-XQ)/dx. 

The  ideal  differential  filter,  F((o),  is  composed 
of  two  filters,  a  linear  amplitude  filter,  A((o)  ,  and  a  phase 
filter,  0 (w) ,  which  are  shown  in  Figure  2.  Since  the 
difference  between  the  negative  and  positive  frequencies  is 
TT  radians  for  the  phase  filter,  it  is  fabricated  as  a  filter 
that  shifts  half  the  frequency  plane  by  tt  radians  [5]  .  The 
linear  amplitude  filter  is  made  using  film,  but  there  is  a 
discontinuity  at  the  origin  since  the  film  will  not  record 

O 

infinite  densities  [4] .  The  anplitude  transmission  of  the 

film  is  limited  to  the  values  A  ,  which  is  less  than 

max 

unity,  and  A^^^^  ,  which  is  greater  than  zero.  The  ampli¬ 
tude  transmission  limits  of  the  film  make  the  linear 
amplitude  filter  bandpass  limited.  The  phase  filter  is 
also  bandpass  limited.  Therefore,  an  approximation  is  made 
of  the  differential  filter;  it  is  the  filter  F' (w)  as  shown 
in  Figure  3  [5] .  The  function  rect  (w/2B)  is  a  high  pass 
filter.  Its  purpose  is  to  eliminate  the  low  spatial  fre¬ 
quencies,  (O' |co|  ^  B)  .  The  amplitude  transmission  value 

of  A  .  occurs  at  the  spatial  frequencies  +  B,  and  the  value 
min  — 

A  occurs  at  the  frequencies  +  W.  The  linear  amplitude 
max  ^  ^ 

filter  has  a  slope,  k,  which  is  generally  less  than  one. 

VJhen  the  input  g(x^)  is  filtered  by  the  approxi¬ 
mate  differential  filter,  the  amplitude,  U(w),  in  the 
frequency  plane  is 
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Figure  3.  The  approximate  differential  filter. 
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U(aj)  =  G(U))F'  (to) 


to  - 


W-B 


(0  + 


W+B 


=  jkio  G((o)  [rect  I — — ]  +  rect  [ — — ]]  (10) 


W-B 


and  the  amplitude  in  the  image  plame  is 


U(Xj^)  =  SF  {G(to)F'  (to)  } 


to  - 


W-B 


=  {jkto  G((o)  }  *  0  {rect  [ — — ] } 


to  + 


W+B 


+  0  {jkto  G(to)}  *  0{rect  [ — — ]} 


(11) 


dg(x.)  - j  (^y ")  X. 

“  k  "d— --  *  (W-B)  sine  [(W-B)x^]e 


dg(x^)  +j(W^)x^ 

+  k  _  ■  (W-B)  sine  [(W-B)x^]e  ^  ^  (12) 


dg(x.)  W+B 

2(W-B)k  ^  ■  *  sine  [  (W-B)  x^]  cos  [  (^)  x^]  (13) 


The  symbol  *  denotes  the  convolution  integral.  This  is  the 
image  amplitude  of  the  input  g(x^)  filtered  by  the  approxi¬ 
mate  differential  filter  F'  (to)  in  a  single  lens  system  such 
as  Figure  1,  page  2. 

The  slope  and  the  bandwidth  of  the  filter  F*  (to)  is 
changed  by  the  magnification  of  the  frequency  plane.  The 
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slope  of  the  filter  F' (w)  is  related  to  the  magnification, 
m,  as 

F'  (w)  =  jmko)  (14) 

and  the  cutoff  frequencies  as 


W 


m 


W 

m 


(15) 


B 


m 


B 

m 


(16) 


The  approximate  one-dimensional  differential  filter  has  a 
complex  linear  amplitude  transmission  along  the  x  axis  and 
constant  amplitude  transmission  along  different  points  on 
the  y  axis.  The  y  axis  is  bandlimited  by  the  function 
rect  ((^/Wy)  just  as  the  x  axis  is  bandlimited  by  the 
function  rect  •  This  does  not  affect  the  derivative 

of  a  one-dimensional  input,  but  it  does  affect  a  two- 
dimensional  input.  As  the  filter  F' (u)  is  rotated  in  the 
frequency  plane,  the  slope  is  decreased  and  the  bandwidth 
changes.  The  filter  F' (u)  in  terms  of  polar  coordinates 
becomes 


F'  (o)  )  =  jk  w 

where 


(17) 


=  k  cos 9 


(18) 
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and  oj^  is  the  spatial  frequency  with  respect  to  the  radius . 
For  0-0° 


F’  (u)^)  =  F*  (w^)  (19) 

The  filter  F' (u^)  will  equal  zero  at  the  angle 


0 


■!  W 
tan-1 

X 


(20) 


With  the  rotation  of  the  filter,  the  upper  cutoff  frequency 
is 


2  .  „  2 


1/2 


Wr  =  (W^"-  +  Wy‘)  cos(0  -  a 


(21) 


where 


-1  w, 

C  =  tan 

X 


(22) 


The  lower  cutoff  frequency  is 


=  B  sec  0 


(23) 


Higher  order  derivatives  are  generated  as  long  as 


a  filter  (jo))'^  is  realized.  The  nth  order  derivative  is  [51 


=  {(ju))^  G(co)} 
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II.  FABRICATION  OF  THE  DIFFERENTIAL  FILTER 

The  Linear  Amplitude  Filter 

The  differential  filter  is  composed  of  two  filters , 
a  linear  amplitude  filter,  A(a)),  and  a  phase  filter,  6  (u)  . 

The  linear  amplitude  filter  is  approximated  by  the  filter 
A' (u)  as  shown  in  Figure  4.  The  filter  A' (w)  is  obtained  on 
film  by  varying  the  exposure  to  the  film.  Since  a  one¬ 
dimensional  amplitude  filter  is  required,  the  exposure  is 
varied  one—iimersionally  by  an  edge  of  a  plate  moving  across 
a  uniform  light  source  at  a  programmed  rate. 

The  filter  A'  (oj)  is  related  to  distance  by 


(25) 


where  x^  is  a  distance  in  the  frequency  plane,  X  is  the 
wavelength  of  the  coherent  light  source,  and  f  is  the  focal 
length  of  a  lens  [1] . 

The  filter  A* (u)  is  a  distance  2d  in  width,  and  d 
is  defined  as 


d  =  XfW  (26) 

Since  the  filter  A' (w)  is  an  even  function,  the  exposure  is 
varied  identically  for  the  distances  +  d.  The  amplitude  of 
the  filter  A' (w)  is  the  amplitude  transmission  of  the  film, 
T^  ,  which  is  a  function  of  exposure.  Hence,  for  the 
correct  variations  of  the  exposure,  the  distance  must  be 
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expressed  as  a  function  of  exposure  to  give  the  correct 
amplitude  transmission  for  the  filter.  The  velocity  of  the 
edge  determines  the  rate-of-change  of  the  exposure,  and  this 
is  derived  from  the  derivative  of  the  distance-exposure 
equation. 


The  correct  exposure  control.  The  amplitude 
transmission  of  the  film  is  a  function  of  the  film  density, 
D,  and  it  is  defined  as  [1] 


The  density  of  the  film  is  a  function  of  the  exposure.  The 
relationship  of  the  exposure  and  film  density  is  best 
expressed  linearly  with  the  well-known  Hurter-Drif field 
curve,  but  this  curve  usually  is  not  linear  for  densities 
less  than  0.8.  Therefore,  the  amplitude  transmission  versus 
exposure  curve  is  used  to  give  a  linear  relationship  of  the 
exposure  and  amplitude  transmission  for  densities  less  than 
0.8.  This  curve  is  linear  for  the  non-linear  portion  of  the 
H  and  D  curve.  Each  of  these  curves  are  shown  in  Figure  5 
for  a  typical  film. 

For  the  H  and  D  curve,  the  film  density  is  [1] 

D  =  Y  log  (|-)  (28) 

o 

where  y  is  the  slope  of  linear  portion,  E  is  the  exposure. 
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(a)  H  and  D  curve 


(b)  Amplitude  transmission 
versus  exposure  curve 


Figure  5.  Film  characteristic  curves. 


AEDC-TR-71-137 


AEDC-TR-71-137 


and  Eq  is  the  value  of  the  exposure  for  the  density  equal  to 
zero.  The  log  is  to  the  base  ten. 

From  Equations  (27)  and  (28) ,  the  amplitude  trans¬ 
mission  of  the  film  is 

E 

Ta(E)  =  (g^)  (29) 

The  amplitude  transmission  versus  exposure  curve 
is  described  as  [10] 


T_  (E)  =  T  -  k,E  (30) 

A  O  1 

where  T^  is  the  amplitude  transmission  for  the  exposure 
equal  to  zero,  and  is  the  slope  of  the  linear  portion. 

For  the  bandpass  section  of  the  amplitude  filter, 

A'  (u)  =  k|a)l  (31) 

or  as  a  function  of  distance  in  the  frequency  plane 

A' (Xj)  =  k' |Xj I  (32) 

where  k'  =  k/Xf .  Since  the  filter  A' (w)  is  made  photo¬ 
graphically. 


T;^(Xf)  =  k'  I 


(33) 
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From  Equation  (29) 


y/2 


K,|  =  (g£) 


(34) 


for  the  linear  portion  of  the  H  and  D  curve,  and  from 
Equation  (30) 


(35) 


for  the  linear  portion  of  the  amplitude  versus  exposure 
curve . 

Both  expressions  for  the  distance,  x^  ,  are 
dependent  on  the  variable  E,  the  exposure.  The  exposure  is 
defined  as  intensity,  I,  times  time,  t(l),  or 

E  =  It  (36) 


This  gives  two  independent  variables.  From  observation,  it 
is  readily  noticeable  that  a  single  independent  variable  is 
much  easier  to  handle,  and  since  time  is  easier  to  vary,  the 
intensity  is  held  constant.  Then  Equation  (34)  becomes 


t-Y/2 


(37) 


and  Equation  (35)  becomes 


'f2 


(38) 
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where  is  a  constant  intensity  and  the  section  of  the 
filter  A'  (w)  based  on  the  H  and  D  curve  is  denoted  by 
and  x^2  denotes  the  section  based  on  the  amplitude  trans¬ 
mission  versus  exposure  curve. 

The  velocity  that  the  edges  must  travel  across  the 
uniform  light  source  to  give  the  correct  exposure  variation 
is  found  by  taking  the  derivative  of  Equations  (37)  and  (38) 
with  respect  to  time.  The  velocity  from  Equation  (37)  is 


V 


1 


^-1/2 (y+2) 


and  the  velocity  from  Equation  (38)  is 


(39)  ■ 


(40) 


The  function  is  valid  for  the  small  values 

of  the  amplitude  transmission  and  the  function  1x^2 1  is 
valid  for  the  larger  values  of  the  amplitude  transmission. 
The  velocity,  v^  ,  is  where 


V 


1 


=  V 


(41) 


The  function,  x^^^  ,  is  invalid  for  ^  ^  snd  the 

function,  x^2  •  i®  invalid  for  V2 

Fabricating  the  linear  amplitude  filter.  The 
method  used  to  traverse  the  edges  must  be  kept  in  mind  when 
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selecting  the  film  to  be  used.  One  of  the  simplest  methods 

is  the  use  of  a  D.  C.  motor  to  move  the  edges.  The  D.  C. 

motor  must  traverse  the  edges  at  the  velocities  required  by 

Equations  (39)  and  (40) .  The  edges  are  mounted  on  threaded 
/ 

rods  with  16  threads/inch.  For  every  16  turns  of  the  rod, 
the  edges  move  one  inch.  The  threaded  rods  are  connected  to 
gears  which  are  linked  with  a  gear  mounted  to  the  motor's 
shaft.  The  ratio  of  the  gears  reduces  the  motor  by  12/19. 
The  method  for  traversing  the  edges  is  discussed  later  in 
this  chapter. 

The  velocity  of  the  edge  is  a  function  of  y  and 

of  the  film  and  of  the  exposure  time.  Gamma,  y,  is  the 

slope  of  the  H  and  D  curve,  k^  is  the  slope  of  the  amplitude 
transmission  versus  exposure  curve,  and  the  exposure  time  is 
the  time  required  to  achieve  a  certain  density.  A  film  with 
a  low  Y  and  kj^  ,  a  high  density  capability,  and  a  long 
exposure  time  is  desired.  The  higher  y  and  k^  ,  the  greater 
the  motor  velocity  range  capability  of  the  motor  needed  to 
move  the  edges.  If  the  film  has  a  low  y,  it  also  has  a  low 
density  capability,  and  vice  versa.  The  shorter  the 
exposure  time,  the  greater  the  velocity  that  the  edges  must 
travel.  Kodak  649-F  film  came  the  closest  to  these  desired 
requirements . 

To  accurately  generate  the  filter  A' (w) ,  the 
characteristic  curves  of  the  film  being  used  must  be  known. 
The  curves  for  the  Kodak  649-F  4x5  inch  plate  was  obtained 
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by  exposing  three  density  step  wedges,  placed  against  the 
film,  at  an  intensity  of  5  foot  candles  for  a  period  of  20 
seconds  and  processing  the  film  in  Kodak  D-19  developer  for 
three  minutes  at  76 “F.  The  resulting  H  and  D  curve  and 
amplitude  transmission  versus  exposure  curve  are  shown  in 
Figure  6. 

From  Figure  6,  the  constants  for  Equations  (39) 
and  (40)  were  calculated  to  give  the  following: 

E^  =  10.0  foot  candles/sec 

y  =  4.57 

T  =  1.22 
o 

=  0.0531 

I  =3.75  foot  candles 
o 

The  velocity  ^2  is  then 


cm/sec 


(42) 


Since  the  velocity  V2  is  the  maximum  velocity  of 
the  edges,  the  slope,  k',  of  the  filter  A' (x^)  is  selected 
to  have  the  velocity  V2  nearest  the  maximvim  capability  of 
the  D.  C.  motor  used  to  traverse  the  edges.  The  slope,  k', 
was  selected  to  be  0.2.  Then,  the  value  of  the  V2  and  v^^  is 


V2  =  +1.0  cm/sec 


(43) 
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=  +  111.5  cm/sec  (44) 

A  plot  of  the  velocities  and  distances  as  a  function  of 
exposure  is  shown  in  Figure  7.  The  velocity  v^  is  easily 
seen  to  occur  at  t  =  4.3  seconds  or  at  an  exposure  of 
16.15  fcs.  Since  the  maximum  density  of  four  could  be 
measured,  the  maximum  exposure  time  is  20  seconds. 

The  edges  travel  from  +0.05  cm  to  +  5  cm  in  a 
period  of  20  seconds,  which  is  shown  in  Figure  7(a).  For 
making  the  filter  A' (x^) ,  the  entire  distance  of  10  cm  is 
imaged  on  the  film  by  a  camera  with  a  magnification  of 
1/lOth.  The  filter 


A'  (Xj)  =  1^1  (45) 

has  a  density  curve  as  shown  in  Figure  8.  The  density  of 
four  from  -0.05  mm  to  +0.05  mm  is  used  to  block  the  spatial 
frequency  of  0  ^  |(i)|  £  B. 

The  velocities  that  the  edges  must  have  to  vary 
the  exposure  are  known  from 'the  velocity-exposure  curve  of 
Figure  7 (b) .  The  apparatus  used  to  move  the  edges  at  the 
correct  velocities  is  called  a  filter  generator  and  it  is 
shown  in  Figure  9.  A  permanent  magnetic  D.  C.  motor  was 
chosen  for  the  filter  generator  since  it  has  a  linear 
relationship  between  the  voltage  input  and  the  motor  speed 
(RPM) .  In  order  to  have  the  correct  velocity-exposure  curve 


26 


(a)  Distance-exposure  curve 


t,  sec 

(b)  Velocity-exposure  curve 


Figure  7.  The  exposure  curves  for  the  filter  A' (x,) . 


Figure  8.  Density  plot  of  the  filter  A' (x^) , 
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Figure  9,  Filter  generator. 
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output  at  the  edges,  a  voltage-exposure  relation  was  found 
in  terms  of  the  input  to  the  D.  C.  motor.  The  input  voltage 
versus  the  RPM  output  of  the  D.  C.  motor  is  shown  in  Figure 
10.  It  is  described  by  the  equation 

=  (0.046) (RPJ)  +  1.4  vdc  (46) 

The  rpms  of  the  D.  C.  motor  are  converted  by  the  filter 
generator  to  the  desired  edge  velocities. 

The  output  of  the  D.  C.  motor  is  geared  down  by 
12/19.  The  edges  are  mounted  on  threaded  rods  with  16 
threads/inch,  and  the  threaded  rods  are  connected  to  the 
gears.  These  threaded  rods  convert  the  rotation  of  the 
motor's  shaft  to  lateral  displacement  of  the  edges.  The 
edges  are  displaced  one  inch  for  every  16  revolutions  of  the 
thread  rods .  The  velocity  of  the  edges  is  expressed  in 
terms  of  the  D.  C.  motor  speed  as 

„  _  tIZt  r(2.54  cm/in.) (RPM) , 
sec)  (16/in. )  ^ 

==  1.675  X  lO”^  RPM  cm/sec 
or 


RPM  =  (5.975  X  10^)v  (47) 

Therefore,  for  the  voltage  input  to  the  motor  in  terms  of 
the  edge  velocity, 
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=  (27.55  V  +  1.4)  vdc  (48) 

From  Equation  (48)  and  Figure  7,  page  26,  it  is  noticed  that 
the  edge  velocities  less  than  0.11  cm/sec  are  beyond  the 
continuous  capabilities  of  the  D.  C.  motor  since  the  break¬ 
away  rpm  of  the  motor  is  56  rpms.  The  edge  velocities  less 
than  0.11  cm/sec  are  obtained  by  pulsing  the  D.  C.  motor 
with  a  voltage  pulse.  The  pulse  motor  control  is  simply  a 
rotating  switch  which  has  four  contacts,  and  when  the  switch 
is  closed,  a  voltage  is  momentarily  placed  across  the  motor 
terminals.  Varying  the  voltage  across  the  switch  will  vary 
the  speed  of  the  motor.  The  curve  for  the  relationship 

a 

between  the  pulse  voltage  and  edge  velocity  is  shown  in 
Figure  11  for  the  switch  contacts  rotating  at  500  rpms. 

The  input  to  the  D.  C.  motor  is  amplified  by  an 

t 

amplifier  with  a  voltage  gain  of  six  and  a  current  gain  of 

4 

approximately  10  .  From  Equation  (48)  and  Figures  7  and  11, 
the  voltage-exposure  curve  is  found,  and  it  is  shown  in 
Figure  12.  The  volcage-exposure  curve  generates  the  edge 
velocity-exposure  curve  of  Figure  7 (b) .  The  input  to  the 
amplifier  is  the  voltage-exposure  curve.  The  period  of  the 
voltage -exposure  curve  is  20  seconds.  This  waveform  is 
generated  by  a  device  called  a  function  generator  which  is 
shown  in  Figure  13.  The  waveform  is  generated  by  placing  a 
voltage  across  a  5000  ohms,  200  watt  variable  wirewound 
resistor.  The  resistor  is  moved  against  a  wire  in  the  shape 
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2,0  4.0  6,0  8,0  10.0  12.H)  14.0  16.0  18.0  20.0 

t,  sec 

Figure  12.  Voltage-exposure  curve  for  fabrication  of  the 
filter  A' (x^) . 
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of  the  voltage-exposure  curve.  Note  presence  of  the  wire  in 
Figure  13.  The  resistor  is  calibrated  to  give  0.25  volts/cm. 
The  resistor  is  moved  20  cm  in  20  seconds  by  a  D.  C.  motor 
and  a  threaded  rod.  The  amplifier,  D.  C.  power  supply 
number  1  which  drives  the  motor  of  the  function  generator, 
and  D.  C.  power  supply  number  2  which  biases  the  resistor, 
compose  the  power  control  unit.  The  block  diagram  demon¬ 
strating  the  manner  in  which  the  filter  A' (x^)  is  generated 
is  shown  in  Figure  14 . 

The  amplitude  transmission  of  the  generated  filter 
A' (x^)  is  found  by  using  a  microdensitometer.  For  the  first 
filter  attempted,  the  density  was  too  high  in  the  center  and 
too  low  at  the  ends  because  of  the  reciprocity  effect.  To 
correct  these  errors ,  the  voltage-exposure  curve  was  changed 
to  the  dashed  curve  shown  in  Figure  12.  Then,  the  filter 
A' (x^)  was  generated  with  the  accuracy  shown  in  Figure  15. 

The  deviation  at  the  ends  is  again  due  to  the  reciprocity 
effect  which  could  not  be  overcome. 

The  Phase  Filter  6  ((jj) 

The  phase  filter  0(w)  is  made  by  depositing  a 
dielectric  coating  on  glass.  The  dielectric  coating  thick¬ 
ness  is  varied  to  give  the  desired  phase  shift.  The  index 
of  refraction  of  air  is  denoted  as  n  ,  and  n.  corresponds 
to  the  dielectric  coating.  The  index  of  refraction  is 
simply  the  ratio  of  the  velocity  of  light  in  a  vacuum,  c,  to 
the  velocity  of  light  in  the  medium,  v,  or  [10] 
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Figure  15.  Amplitude  transmission  versus  distance  for  the 
experimental  filter  A' (x^) . 
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(49) 


As  seen  from  Figure  2,  page  10,  a  phase  of  -  Tr/2  radians  and 
+  Tr/2  radians  is  needed  for  the  fabrication  of  the  filter 
0(uj)  which  is  a  total  phase  difference  of  tt  radians  between 
the  positive  and  the  negative  half  frequency  plane.  There¬ 
fore,  a  phase  shift  of  tt  radians  is  only  needed  for  the 
negative  frequency.  For  the  phase  shift  of  it  radians,  the 
light  wave  must  be  delayed  in  the  dielectric  material  by  X/2 
with  respect  to  light  wave  traveling  in  air.  Hence, 


Ub)f 


b  a 

where  t'  is  the  time  required  for  the  light  to  travel  a 
distance  X/2  in  the  dielectric  material.  Then  the  required 
thickness  of  the  dielectric  material  is 


d 


n 


(51) 


Fabricating  the  phase  filter.  The  required  thick¬ 
ness  of  the  dielectric  coating  is  deposited  on  one  half  of  a 
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one-inch  square  glass  plate.  To  measure  the  amount  of 
dielectric  deposited  on  the  glass,  the  reflectivity  of  the 
dielectric  material  on  the  gl^ss  is  measured.  The  reflec¬ 
tivity  as  a  function  of  wavelength  thickness  will  resemble 
Figure  16  [10]  .  The  one-inch  square  glass  plate  is  placpci 
on  top  of  a  stand  inside  a  bell-shaped  vacuum  chamber.  One- 
half  of  the  glass  plate  is  blocked  by  a  sheet  of  metal.  The 
dielectric  material  is  deposited  by  heating  it  to  a  high 
temperature  at  the  bottom  of  the  stand.  Due  to  the  shape  of 
the  vacuum  chamber,  the  light  source  used  to  measure  the 
thickness  of  the  dielectric  is  incident  upon  the  glass  plate 
at  45®.  The  light  source  is  a  He-Ne  laser.  Therefore,  the 
thickness  measured  by  the  reflectivity  is  at  an  angle  of  71/4 
radians.  Hence,  the  normal  thickness,  d^^  ,  is 

d  =  0.707  d  (52) 

n  it/4 

The  dielectric  material,  silicon  mono-oxide,  SiO,  is  used. 
Its  index  of  refraction  is 


%  =  1.97 


then 


(53) 


d^  =  0.516\ 


(54) 
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Reflectivity 
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Since  only  thickness  is  measured,  the  dielectric 

coating  is  stopped  when  the  thickness 


d^/4  =  0.74X 


(55) 


is  reached. 


42 


AEDC-TR-71-137 


CHAPTER  III 

EXPERIMENTAL  COMPARISONS  OF  VARIOUS 
OPTICAL  FILTERING  METHODS 

I .  INTRODUCTION 

Experiments  are  constructed  to  test  the  charac¬ 
teristics  of  various  optical  filters.  The  diagram  of  the 
test  apparatus  is  shown  in  Figure  17  (a) .  The  TV  system  is 
used  to  view  the  image  plane;  it  is  a  high  resolution  system 
of  a  1000  horizontal  scan  lines  per  vertical  scan  with  a 
bandwidth  of  32  MHz.  To  examine  closely  the  effect  of  • 
optical  filtering  on  the  input,  a  portion  of  a  single  hori¬ 
zontal  scan  line  of  the  TV  system  is  viewed  on  an 
oscilloscope;  this  relates  to  a  cross  section  of  the 
intensity  pattern  in  the  image  plane.  The  oscilloscope  used 
is  a  Tektronix  547  with  a  1A5  preamp  yielding  an  oscillo¬ 
scope  rise  time  of  8  nsec.  Therefore,  the  combination  of 
the  oscilloscope  and  the  TV  system  has  a  rise  time  of  13.6 
nsec.  The  TV  system  was  calibrated  so  that  the  distance 
could  be  measured  from  the  time  scale  of  the  oscilloscope. 
The  relationship  is  131.0  +4.65  microns/cm  for  a  time 
scale  of  0.2  ysec/cm  for  unity  magnification. 

The  ability  of  the  various  filters  to  minimize  the 
error  in  the  location  of  an  edge  is  tested  using  a  step 
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(a)  Diagram 


Figure  17.  Filter  test  apparatus 
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(b)  Photograph  of  experiment 


Figure  17.  (continued) 
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function.  In  addition,  several  of  the  filters  are  applied 
to  phase  information  for  comparison  of  the  various  filtering 
techniques . 


II.  THE  STEP  FUNCTION 

The  input  used  to  test  the  response  of  the 
different  optical  filters  is  a  one-dimensional  step  function 
which  is  the  edge  of  a  razor  blade.  The  step  function  has 
an  amplitude  transmission  of 


Then, 


s  (Xo)  =  j  (1  +  sgn  x^) 


(56) 


GF  {s  (x^) }  =  S  (io) 


=  I  [6(0))  +  (57) 

The  step  function  is  imaged  on  the  vidicon,  but  it  is  not  a 
perfect  step  function  due  to  bandlimiting  by  the  optics. 

The  image  intensity  of  the  step  function  is  [12] 

ll(x)  =  l|+  isi(^)l  (58) 

where  q  is  the  image  distance,  b  is  the  radius  of  the  lens, 
and  K  =  2it/X.  The  curve  generated  by  Equation  (58)  is  shown 
in  Figure  18. 
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Intensity 


Figure  18.  Image  intensity  of  a  step  function  bandlimited 
by  a  lens . 
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A  portion  of  a  horizontal  scan  line  showing  the 
intensity  cross  section  of  the  step  function  in  the  image 
plane  is  shown  in  Figure  19.  The  rise  time  of  the  step 
function  is  80  nsec  which  is  slower  than  the  rise  time 
capability  of  the  electronic  system.  This  indicates  that 
the  electronic  system  gives  a  reasonably  accurate  display  of 
the  signals  in  the  image  plane.  The  rise  time  of  the  step 
function  is  therefore  due  to  bandlimiting  effects  as 
described  by  Equation  (58) . 

For  coii^arison  of  the  different  optical  filters, 
the  manual  gain  of  the  TV  system  was  maintained  constant. 
Referring  to  Figure  17(a),  page  44,  lens  L2  magnifies  the 
transform  plane  of  lens  LI  by  0.845,  but  the  overall  magni¬ 
fication  of  the  object  plane  by  lens  Ll  and  L2  is  1.79.  A 
He-Ne  laser  is  used  for  the  coherent  light  source;  therefore, 
in  the  frequency  plane  of  lens  Ll,  each  centimeter  repre- 

3 

sents  a  spatial  frequency  of  4.13  x  10  lines/mm  and  in  the 

3 

frequency  plane  of  lens  L2,  3.49  x  10  lines/mm.  Due  to  the 
magnification  of  the  object  plane,  the  relationship  of  the 
distance  and  time  base  is  33.8  +1.3  microns/cm  for  an 
oscilloscope  time  scale  of  0.1  ysec/cm,  respectively. 

The  step  function  is  attenuated  by  a  neutral 
density  filter  with  an  intensity  transmission  of  0.0467; 
hence,  the  amplitude  of  the  step  function  is  4.28v  instead 
of  200  mv  as  shown  in  Figure  19(a).  These  values  of 
voltages  are  read  off  the  oscilloscope  and  are  proportional 
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y  axis:  100  mv/cm  x  axis:  0.1  p,sec/cm 

(a)  A  horizontal  scan  line  of  the  image  plane 
attenuated  by  a  neutral  density  filter 


(b)  The  image  plane  as  viewed  on  the  TV  monitor 
Figure  19.  The  bandlimited  step  function. 
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I 


to  the  intensity  incident  upon  the  vidicon  active  surface. 
They  are  relative  values  only. 

III.  THE  DERIVATIVE  OF  THE  STEP  FUNCTION 

The  derivative  of  the  step  f motion  is  taken  by 
placing  the  phase  filter,  0(w),  in  the  frequency  plane  of 
lens  Ll  and  the  linear  amplitude  filter,  A'  (to)  ,  in  the 
frequency  plane  of  lens  L2  to  yield  the  filter  F' (u) .  The 
filter  A' (u)  has  dimensions  of  7  mm  x  7  mm,  and  the  phase 
filter  0  (ti))  has  dimensions  of  2.4  cm  x  2.4  cm.  The  spectr\im 
of  the  frequency  plane  of  lens  L2  is  shown  in  Figure  20. 

From  Equation  (13)  ,  the  intensity  of  the  bandlimited 
derivative  of  a  step  function  is 

2 

KxJ  =  4  A^^  k^(W-B)^  [sine  ((W-B)x.) cos  ((^^)x.)l  (59) 

1  max  X 

A  portion  of  the  horizontal  scan  line  in  Figure 
21(a)  shows  the  intensity  of  the  bandlimited  derivative;  its 
maximum  amplitude  is  130  rav.  The  rise  time  of  the  deriva¬ 
tive  is  slower  than  the  rise  time  of  the  step  function. 

This  is  due  to  the  derivative  having  an  upper . bandlimi ting 
frequency,  W.  This  frequency  is  much  smaller  than  the 
limiting  frequency,  b,  of  the  step  function.  The  band- 
limited  derivative  of  the  step  function  is  108  microns  wide. 
The  maximum  error  in  locating  the  edge  is  54  microns. 

The  manual  gain  of  the  TV  system  was  increased, 
and  the  derivative  of  the  step  function  is  shown  in 
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Figure  20.  The  frequency  spe 
filtered  by  the  filter  F' 
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y  axis:  20  mv/cm  x  axis:  0.1  p.sec/cm 

(a)  A  horizontal  scan  line  of  the  image  plane 


y  axis:  100  mv/cm  x  axis:  0.1  (j.sec/cm 

(b)  A  horizontal  scan  line  of  the  image  plane  for 
the  manual  gain  increased 

Figure  21.  The  bandlimited  derivative  of  the  step  function. 
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(c)  The  image  plane  as  viewed  on  the  TV  monitor 


Figure  21.  (continued) 
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Figure  21(b).  The  increase  in  gain  increases  the  amplitude 
of  the  signal  and  increases  the  signal-to-noise  ratio. 

The  filter  F' (u)  is  rotated  in  the  frequency  plane 
and  the  derivative  of  the  step  function  is  shown  to  conform 
to  Equation  (17)  by  Figure  22.  It  is  impossible  to  demon¬ 
strate  the  derivative  as  a  function  of  the  magnification 
since  the  amount  of  energy  per  unit  area  changes  in  the 
image  plane  as  a  function  of  the  magnification  which  causes 
changes  in  the  signal  amplitude. 

IV.  FILTERING  OF  THE  STEP  FUNCTION  BY  THE 
LINEAR  AMPLITUDE  FILTER 


The  filter  A*  (w)  was  discussed  in  Chapter  II.  Its 
amplitude  transmission  is  shown  in  Figure  4,  page  16.  The 
filter  A'  (w)  filters  the  step  function  in  the  frequency 
plane  of  lens  L2;  the  frequency  spectrum  is  shown  in  Figure 
23.  This  frequency  spectrum  is  simply 


0)  + 


W+B 


U(w)  =  k  A  rect  [- 
max 


W-B 


,  ■  IT 

2-, 

]e 


0)  - 


W-B 


+  k  A  rect  ( — r,  _ 
max  W-B 


■  JL 

2  ,  "^2 
le 


(60) 


Then  the  amplitude  in  the  image  plane  is 
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Figure  22.  The  bandlimited  derivative  of  the  step  function 
as  a  function  of  the  angle  0. 
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U(ai) 


i 


Figure  23.  The  frequency  spectrum  for  the  step  function 
filtered  by  the  filter  Cwl . 
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3  {U(t0)  }  =  U(X^) 


“D 

=  3k  (W-B)  [since  (W-B)^Je 

max  XI 

•  /W+B. 

+D  (-Y-) 

-  sine [ (W-B) x^]e  ] 

=  2A  k  (W-B)  [since  (W-B)  X.  1  sin [  (5^^)  X.  ]  ]  (61) 

max  1  z  1. 

The  intensity  in  the  image  plane  is 

2 

I(x^)  =  4k^  A^^^(W-B)  ^  [since  (W+B)d^] sin [  1  (62) 

Notice  that  the  only  difference  between  Equations  (62)  and 
(59)  is  a  sine  and  cosine  fvinction. 

The  intensity  of  the  step  function  filtered  by 
A' (w)  is  seen  in  Figure  24.  The  "valley"  in  the  horizontal 
scan  line  is  caused  by  the  sine  function  in  Equation  (62) 
since  the  sine  is  zero  at  =  0.  This  "valley"  is  used  to 
locate  the  edge  of  the  razor  blade.  The  "valley"  is  40.7 
microns  wide;  hence,  the  error  in  locating  the  edge  is  20.35 
microns.  The  maximum  amplitude  of  Equation  (62)  in  Figure 
24  is  150  mv. 
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y  axis:  50  mv/cm  x  axis:  0.1  p.sec/cm 

(a)  A  horizontal  scan  line  of  the  image  plane 


(b)  The  image  plane  as  viewed  on  the  TV  monitor 
Figure  24.  The  step  function  filtered  by  the  filter  A' (w) . 
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V.  THE  HILBERT  TRANSFORM  METHOD  APPLIED 
TO  THE  STEP  FUNCTION 

The  Hilbert  transform  is  defined  as  [2] 


,  f  U^(x') 

U  (X  )  =  i  I  -7-2 — rr  dx' 

O  O  TT  J  (X^-X') 


(63 


The  Hilbert  transform  filter  is  described  by 


H(a))  =  -j  sgn(w) 


(6) 


This  filter  is  constructed  the  same  as  the  phase  filter, 

6  (u)  .  Therefore,  the  filter,  6(10),  is  used  to  take  the 
Hilbert  transform  of  a  function  in  the  object  plane.  For 
finding  the  Hilbert  transform  of  the  step  function,  the  step 
function  is  represented  by  a  rectangular  pulse  in  the  object 
plane;  this  is  allowed  since  the  collimated  coherent  light 
incident  upon  the  object  plane  is  not  infinite  in  area. 

Then  the  Hilbert  transform  of  the  rect  function 

X 

s'(x^)  =  A  rect  (:p)  (64) 

is 


i'  =  I  In  - ^1  (65) 

Xi  t  ^ 

/\ 

s' (x  )  denotes  the  Hilbert  transform.  The  plot  of 
o 
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Equation  (65)  is  shown  in  Figure  25.  The  edge  of  the  razor 
blade  is  located  at  -t/2.  The  portion  of  the  horizontal 
scan  line  showing  the  edge  at  -t/2  in  the  image  plane  is 
shown  in  Figure  26 (a) .  The  width  of  the  Hilbert  transform 
of  the  razor  blade  edge  is  330  microns.  This  appears 
visually  as  a  smear  of  the  razor  blade  edge.  The  maximum 
error  in  locating  the  edge  is  165  itdcrons.  The  signal 
amplitude  is  attenuated  by  a  neutral  density  filter  with 
intensity  transmission  of  0.0467;  therefore,  the  maximum 
signal  amplitude  is  5.35  volts. 

VI.  THE  SCHLIEREN  METHOD  APPLIED  TO 
THE  STEP  FUNCTION 

Filtering  by  the  Schlieren  method  is  achieved  by 
introduction  of  a  razor  blade  edge  into  the  frequency  plane 
to  block  one  half  the  frequency  spectrum.  This  filter  is 
applied  to  the  step  function  in  the  frequency  plane  of  lens 
L2.  The  Schlieren  method  filter  is  described  as 

T(to)  =  +  sgn{(*)))  (4) 

When  this  filtering  is  applied  to  the  step  function,  the 
image  plane  amplitude  is  [1] 


Ui(Xi)  =  + 


s' (x' ) 

-7-2 - rr  dx’  ] 

(Xj^-X'  ) 
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y  axis:  50  mv/cm  x  axis:  0.1  iisec/cm 

(a)  A  horizontal  scan  line  of  the  image  plane 
attenuated  by  a  neutral  density  filter 


(b)  The  image  plane  as  viewed  on  the  TV  monitor 
Figure  26.  The  Hilbert  transform  of  the  step  function. 
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or 

Ui(Xi)  = 

The  image  intensity  is 

li(Xi)  =  Ttsi^(Xl)  +  (67) 

When  viewing  a  portion  of  a  horizontal  scan  line 
showing  the  edge  of  the  step  function,  the  first  term  of 
Equation  (67)  is  not  detected  since  it  is  a  constant.  Only 
the  second  term  of  Equation  (67)  is  detected  and  it  is  shown 
in  Figure  27(a).  The  maximum  image  intensity  amplitude  is 
500  mv.  The  edge  of  the  step  function  appears  as  an 
intensity  pattern  with  a  width  of  169  microns.  The  maximum 
error  in  locating  the  edge  is  84.5  microns. 

yil.  BANDPASS  FILTERING  OF  THE  STEP  FUNCTION 

.  The  bandpass  filter  is  composed  of  two  filters,  a 
high  frequency  pass  filter  and  a  low  frequency  pass  filter. 
The  high  pass  filter  is  a  number  20  AWG  wire  placed  through 
the  origin  of  the  frequency  plane.  The  low  pass  filter  is 
effected  by  the,  finite  diameter  of ’the  lens.  The  bandpass 
filtering  of  the  step  function  is  performed  in  the  fre¬ 
quency  plane  of  lens  L2.  The  diameter  of  the  lenses  is 
7.6  cm.  The  equation  for  the  amplitude  of  the  bandpass 
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y  axis:  100  mv/cm  x  axis:  0.1  |isec/cin 

(a)  A  horizontal  scan  line  of  the  image  plane 


(b)  The  image  plane  as  viewed  on  the  TV  monitor 

Figure  27.  The  step  function  filtered  by  the  Schlieren 
method . 
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filtered  step  function  is  derived  by  Fujimura  [13] .  The 
equation  is 

G  (x)  =  ^[Si(2iTbx)  -  Si{2iTax)]  (68) 

IT 

For  unity  magnification  where  a  is  the  lower  spatial  cutoff 
frequency  and  b  is  the  upper  spatial  cutoff  frequency.  G(x) 
is  the  amplitude  in  the  image  plane  and  it  is  recognized  as 
an  odd  function.  The  amplitude  of  the  positive  half  plane 
is  plotted  in  Figure  28  for  different  values  of  b.  The 
bandpass  filtered  step  function  appears  visually  as  a 
diffraction  pattern  which  is  symmetrical  about  the  point 
X  =  0.  The  position  of  the  horizontal  scan  line  showing  the 
intensity  of  the  bandpass  filtered  step  function  in  the 
image  plane  is  shown  in  Figure  29.  The  "valley"  in  the 
bandpass  filtered  step  function  occurs  at  x  =  0.  It  is  used 
to  locate  the  edge  of  the  step  function.  The  width  of  the 
"valley"  is  47.4  microns.  The  maximvim  error  in  locating  the 
edge  is  23.7  microns.  The  signals'  maximum  intensity  is 
500  mv. 

Vill.  APPLICATIONS  OF  THE  VARIOUS  FILTERS 

A  summary  of  the  operation  of  the  previously  con¬ 
sidered  filters  upon  a  step  function  and  the  ability  to 
locate  object  edges  is  presented  in  Table  I.  The  various 
optical  filters  were  applied  to  different  applications  to 
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Figure  28.  Amplitude  in  a  bandpass  filtering  image  of 
an  Edge. 
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y  axis:  100  mv/cm  x  axis:  0.1  p.sec/cm 

(a)  A  horizontal  scan  line  of  the  image  plane 


(b)  The  image  plane  as  viewed  on  the  TV  monitor 
Figure  29.  Bandpass  filtering  of  the  step  function. 
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TABLE  I 

SUMMARY  OF  THE  EFFECTS  OF  VARIOUS  FILTERS  ON  THE  STEP  FUNCTION 

0 


Step  Function  and 
Filter 

Maximum  Intensity 
of  Signal  in 
the  Image  Plane 

Percentage  of 
Attenuation 

Maximum  Error 
for 

Locating  Edge 

Unfiltered 

4.28 

4 

volts 

0% 

Differential  Filter, 
F'  (oi) 

0.13 

volt 

97.0% 

54.00 

microns 

Linear  Amplitude 
Filter,  A' (w) 

0.15 

volt 

96.5% 

20.35 

microns 

Hilbert  Transform 

5.35 

volts 

-25.0% 

169.00 

microns 

Schlieren  Method 

0.5 

volt 

88.3% 

84.5 

microns 

Bandpass  Filter 

0.5 

volt 

88.3% 

23.7 

microns 
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provide  an  insight  into  their  operation.  The  first  example 
is  the  wire  mesh  shown  in  Figure  30(a).  The  wire  mesh  is  a 
two-dimensional  object.  In  Figure  30(b)  the  wire  mesh  is 
filtered  by  the  phase  filter  0(io).  The  filter  6  (lo)  is  one¬ 
dimensional.  It  yields  the  Hilbert  transform  of  the  wire 
mesh  only  for  the  x-coordinate .  The  "smear"  effect  is 
easily  observed.  Notice  the  y -coordinates  are  unaffected. 

In  Figure  30(c)  the  wire  mesh  is  filtered  by  the  approximate 
one-dimensional  filter  F' (w) .  Again,  only  the  x-coordinates 
are  affected.  The  derivative  is  taken  of  both  edges  of  the 
vertical  wires. 

The  second  example  is  the  hologram  of  an  air 
stream  escaping  at  a  high  velocity  from  a  nozzle  into  the 
flow  field  of  a  wind  tunnel.  The  disturbances  in  the  center 
of  Figure  31(a)  is  the  air  stream  from  the  nozzle.  This 
unfiltered  air  stream  shows  very  little  phase  information. 

In  Figure  31(b)  the  flow  field  is  filtered  using  the  one¬ 
dimensional  Hilbert  transform  filter.  The  phase  information 
becomes  very  prominent,  especially  the  air  flow  from  the 
nozzle.  Notice  the  large  shock  wave  which  extends  from  the 
top  to  the  bottom  of  the  figure.  This  shock  wave  is  not 
visible  in  Figure  31(a)  since  it  is  a  small  phase  change. 

In  Figure  31  (c)  the  flow  field  is  filtered  by  the  one¬ 
dimensional  Schlieren  method.  The  large  shock  wave  is 
visible,  but  the  phase  information  around  the  nozzle  is  not 
visible.  The  phase  information  is  affected  by  the 
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(a)  Unfiltered 


(b)  Filtered  by  the  filter  6  (o)) 
Figure  30.  'Filtering  of  a  wire  mesh. 
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(c)  Filtered  by  the  filter  F' (w) 
Figure  30.  (continued) 


\ 
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(a)  Unfiltered  (flow  is  from  left  to  right) 


(b)  Filtered  by  the  filter  9 (w) 
Figure  31.  Filtering  of  a  flow  field. 
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(c)  Filtered  by  the  Schlieren  method 


Figure  31.  (continued) 
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one-dimensional  filter  as  long  as  it  has  an  x-coordinate . 
Filtering  of  the  flow  field  by  the  differential  filter  F' (w) 
was  attempted,  but  the  filter  F*  (uj)  greatly ■  attenuated  the 
information  of  the  hologram,  and  no  useful  data  were 
obtained. 

The  third  example  is  an  air  stream  from  a 
compressed  air  gun  striking  the  end  of  a  rod.  For  the 
unfiltered  air  stream,  no  phase  information  is  visible. 

When  the  air  stream  is  filtered  by  the  Hilbert  transform 
filter,  the  phase  information  becomes  visible. as  shown  in 
Figure  32 (b) .  The  air  stream  is  filtered  by  the  differ¬ 
ential  filter  F' (w)  in  Figure  32(c).  The  differential 
filter  accentuates  the  phase  information  that  is  visible  in 
Figure  32  (b) .  The  one -dimensional  differential  filter 
affects  all  information  in  the  x  direction. 
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(a)  Unfiltered  (flow  is  from  right  to  left) 


(b)  Filtered  by  the  filter  0 (u) 

Figure  32.  Filtering  of  a  compressed  air  stream. 
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(c)  Filtered  by  the  filter  F’ (u) 
Figure  32.  (continued) 
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CHAPTER  IV 

CONCLUSIONS 

In  this  study  a  comparison  was  made  between  the 
differential  filter,  F' (w) ,  the  linear  amplitude  filter, 

A' (w)  ,  the  Hilbert  transform  filter,  0(a)),  the  Schlieren 
method,  and  bandpass  filtering.  Each  filter  was-  fabricated 

and  tested.  In  addition",  each  of  the  various  filters  was 

,  1 

applied  to  a  step  function  to  observe  which  filter  was  best 
for  locating  the  edge  of  the  step.  Two  conditions  had  to  be 
taken  into  account  when  considering  the  various  filters,  the 
intensity  of  the  filtered  step  function  and  the  error  in 
locating  the  edge. 

It  was  found  that  the  Hilbert  transform  method,  or 
the  phase  filter  9  (o))  ,  had  too  great  an  error  in  locating 
the  edge,  a  169  micron  error.  The  Schlieren  method  had  an 
error  of  84.5  microns  which  was  also  too  great  when  con^ared 
to  the  bandpass  filter  error  of  23.7  microns.  The  linear 
amplitude  filter.  A' (o))  ,  had  the  minimum  error  of  20.35 
microns  in  locating  the  edge;  however,  its  intensity  was  too 
low. 

The  intensity  level  of  the  step  function  filtered 
by  the  differential  filter,  F' (w) ,  was  the  lowest.  However, 
the  differential  filter  had  an  error  of  54  microns  in 
locating  the  edge  which  was  twice  as  great  as  the  error  of 
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the  bandpass  filtering  method  and  also  of  the  linear  ampli¬ 
tude  filter.  Therefore,  it  was  concluded  that  the  bandpass 
filter  was  the  best  filtering  method  for  locating  the  edge 
of  the  step  function;  in  addition,  except  for  the  phase 
filter,  6  (u) ,  it  has  the  highest  intensity  level  in 
filtering  the  step  function. 

Next,  the  Hilbert  transform  filter,  the  Schlieren 
filter,  and  the  differential  filter  were  compared  in  appli¬ 
cations  to  test  their  ability  to  identify  phase  information. 
The  first  application  of  the  flow  field  in  the  wind  tunnel 
brought  out  the  limitations  of  the  differential  filter  since 
it  could  not  be  used  on  a  low  intensity  signal.  When  a  more 
intense  signal  was  used,  the  differential  filter  produced 
better  results  comparable  to  that  of  the  Hilbert  transform. 
The  Hilbert  transform  filtering  of  the  flow  field  clearly 
identified  the  phase  information.  When  the  flow  field  was 
filtered  by  the  Schlieren  method,  half  the  phase  informa¬ 
tion  was  obscured  by  noise,  producing  very  poor  results. 
Thus,  it  appears  in  some  cases  that  the  best  filter  for 
identifying  phase  information  is  the  Hilbert  transform 
filter. 
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